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Abstract 



We call a state "vacuum-bounded" if every measurement performed outside a 
specified interior region gives the same result as in the vacuum. We compute 
the maximum entropy of a vacuum-bounded state with a given energy for a 
one-dimensional model, with the aid of numerical calculations on a lattice. 
The maximum entropy is larger than it would be for rigid wall boundary con- 
ditions by an amount 5S^ which for large energies is ^ | In(LinT), where Lin 
is the length of the interior region. Assuming that the state resulting from 
the evaporation of a black hole is similar to a vacuum-bounded state, and 
that the similarity between vacuum-bounded and rigid-wall-bounded prob- 
lems extends from 1 to 3 dimensions, we apply these results to the black hole 
information paradox. We conclude that large amounts of information cannot 
be emitted in the final explosion of a black hole. 
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I. INTRODUCTION 



A. Background 

Since the discovery of black hole radiation by Hawking , the fate of information falling 
into a black hole has been a mystery. (See for reviews.) If Hawking's semiclassical 

calculation is correct, then the outgoing radiation is purely thermal and the outgoing photons 
are uncorrelated to each other and to the matter which formed the black hole. If the 
evaporation is complete, and if the thermal nature of the radiation persists throughout 
the evaporation, then the original information is lost. Thus if the black hole is formed 
from a quantum-mechanically pure state, there will nevertheless be a mixed state after the 
evaporation is complete. This is the position held by Hawking (e.g. see |^), but it violates 
CPT and may lead to difficulties with energy conservation and causality 0J^J^,^. 

If information is not lost in black hole evaporation, there are several possibilities. One 
is that the black hole does not evaporate completely, but instead produces one or more 
Planck-scale remnants (e.g. see [^). Another possibility is that the information disappears 
into a baby universe [Q. In this scenario the quantum- mechanical pure state is preserved, 
but parts of it are inaccessible to observation. It is also possible that the radiation is not 
really thermal, even at early times, because of a complementarity principle |llO| , [Tl[| or the 



inapplicability of the semiclassical approach p!^JT3|JT^ , and thus that the information is 
encoded in subtle correlations in the radiation. In this case the black hole could act like a 
normal object with the entropy describing internal degrees of freedom. Some results from 
string theory ||l5| , p!6| tend to confirm this view. 



Even if the radiation is thermal and uncorrelated during most of the evaporation, there 
is no reason to believe that it remains thermal near the endpoint of the evaporation. The 
late-time radiation is presumably governed by an unknown theory of quantum gravity, and 
may well have correlations to the radiation emitted earlier .0 However, it is generally believed 



that late-time radiation cannot resolve the information paradox [p77|,|3|,|4| . The argument goes 
as follows: While the black hole is large, it is presumably radiating high-entropy thermal 
radiation. If the final explosion is to restore a pure state, it must radiate as much entropyll 
as was radiated in earlier times. However, by the time the black hole reaches the point where 
unknown physics could come into play, there is little energy remaining. To radiate a lot of 
information with little energy requires a long period of time, and thus the "final explosion" 
looks more like a long-lived remnant. 



However, Wilczek |TB| argues from a moving mirror model that a state with high entropy 
can nevertheless be purified with arbitrarily low energy cost. In certain ways, his model looks 
more like a remnant theory than a complete-evaporation theory, but it still appears to cast 
some doubt on the standard argument above. 



^This is conceivable because at early times information in the outgoing radiation can be correlated 
with information in the ingoing negative-energy flux. 

^Here and throughout this paper, "entropy" means fine-grained quantum-mechanical entropy. 
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In any case, this argument requires bounding the entropy that can be contained in a 
particular region with a fixed energy. In the case of a region with refiecting walls, this is the 
question of finding the thermal state of quantum fields in a box. For a spherical box and a 
particular field theory the problem is easily solved. But with a region of complex shape, or 
where one wishes to make a statement intended to apply to all field theories, the situation 
is more complicated. Bekenstein argues that such a universal bound exists, although 



Unruh and Wald |2T[] disagree. 

Here we take a different approach. We consider only a single scalar field, but we use a 
weaker and, we hope, more physical condition on the results of the black hole evaporation. 
In the end, our results still support the claim that late-time radiation cannot restore the 
purity of the state of an evaporating black hole. 



B. The "vacuum-bounded state" 

We start by considering a black hole formed from a pure quantum-mechanical state of 
incoming matter. To avoid any possible complications of quantum gravity theory, we will 
look at the state produced after the black hole has completely evaporated 0. Gravity 
should play no significant role in this state, since the energy density state should be small 
everywhere.! We can describe the final state as follows: At large distances from the position 
of the black hole (taken as the origin) there is outgoing Hawking radiation, which we are 
assuming to be thermal and uncorrelated to anything. Within some distance R of the 
origin, there is some state of ordinary quantum fields that could have correlations with the 
radiation emitted earlier. The distance R is the distance to which such information might 
have propagated after unknown physics came into play. Let us assume that Hawking's semi- 
classical calculation is good up to an energy scale Tunk- This temperature is reached when 
the black hole's mass isQ Mq = l/(87rTunk)- If, after this, the rate of evaporation continues 
to match the Hawking calculation,i the black hole will evaporate in time t ~ 10^ M^/g ~ 
l/{gT^^^^), where g is the effective number of degrees of freedom in the particles that can be 
radiated. (See ||2^.) So there is a sphere of radius 

R-z^ (1-1) 

y-^unk 



which contains total energy 



Eo = —— (1.2) 

onl unk 



^If instead there are Planck-scale concentrations of energy, then we would have a remnant theory, 
a possibility we are explicitly not considering here. 

^We are working with units in which c = G = h = kB = l 

^As opposed, for example, to slowing to nothing and leaving a remnant. 



3 



in which the information could be contained. Taking, for example, Tunk = lO^^Gev ~ 
10~^mpi and g ~ 100 we geti 

R ~ 10^°/pi ~ lO-^^cm (1.3a) 
Eo ~ lO^rripi ~ lO^^g ~ lO^^erg . (1.3b) 

Now we would like to answer the following question: How much entropy can be contained 
in a spherical region of radius R with a total energy of E'o? To answer this question we have 
to specify what we mean by "contained in a region". As mentioned earlier, if we ask how 
much entropy can be contained in a spherical box of radius R with perfectly reflecting walls, 
the question can be easily answered. However, the system with the box is not so closely akin 
to the system under discussion. For instance, inserting the reflecting walls into the system 
produces a divergent increase in the ground-state energy of the system. Furthermore, if 
we started with the vacuum in the whole system, and then introduced a spherical wall, we 



would produce a divergent geometric entropy p3| , p^ . A better description of our system is 



simply that it has thermal radiation outside radius R, and an unknown state of the quantum 
fields inside radius R, but no barrier or boundary at R. 

To study such systems we will assume that the difference between the external Hawking 
radiation and an external vacuum is not important to considerations of entropy.i We will 
study systems that have an arbitrary state inside R but the vacuum outside R. To make 
this precise we will specify the problem as follows: 

Let a "vacuum-bounded state" be a generalized state (i.e. density matrix) 
for which every operator composed of fields at points outside a specified interior 
region has the same expectation value as in the vacuum. What is the maximum 
entropy of such a state whose interior region is a sphere of radius R and whose 
total average energy§ is given by (H) = EqI 

We expect to find that the answer to this question is similar to that of a box of radius 
R with reflecting walls, with some small correction. 



II. GENERAL CONSIDERATIONS 

We will first consider a general system divided into an inside region ("part 1") and 
an outside region ("part 2"). We will say that a generalized state (i.e. density matrix) 



^Another possibility is that g diverges as T ^ mpi. In this case the information can be radiated 
in a small number of particles of about the Planck mass, chosen from an infinite spectrum of such 
particles. This is effectively a remnant theory. 

^If this approximation is bad we can always increase R until the Hawking radiation outside R has 
very low temperature. 

*We cannot specify that any measurement of the energy must give a particular value. Such a 
state is necessarily static and thus cannot represent outgoing radiation. 
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is "localized to the inside" or "obeys the vacuum-bounded condition" if any measurement 
performed on the outside field operators in this state yields the same result as in the vacuum, 
i.e. if 

TrpOs = Trp™^02 = (OlOslO) (2.1) 

for every operator O2 that is constructed out of field operators in the outside region. 

In the language of density matrices, we can write p2 = Tri p, where p is the overall 
density matrix describing our system and Tri means to trace over all the "inside" variables. 
Then p2 is the reduced density matrix describing only the outside variables. Eq. (|2.1| ) then 
is equivalent to 

P2=pr = Tri|0)(0|, (2.2) 
where |0) denotes the ground state. 



A. Stationary points of S 

The present problem is a particular case of the following general problem: Find a density 
matrix p (i.e. a positive-semidefinite Hermitian matrix with unit trace) which maximizes 
S = — Trplnp subject to a set of constraints Tr pC^ = V^. To analyze this problem we will 
look at the change in S{p') when p is varied via p' = p + t6p. The naive result holds for the 
first derivative, 

— = -Tr ((5p In p + (5p) . (2.3) 

(Jjv 

We will only consider variations of p which leave p normalized with Tr p = 1 , which means 
that Tr5p = and thus 

— = -Tropin p. (2.4) 
at 



We show in appendix A that entropy is always strictly (downward) concave ||2^ which 
means that 

for any 6p. Thus any point at which S is stationary is a local maximum. Furthermore, we 
can show that there can be at most one such maximum, as follows. 

Let V be the space of all Hermitian matrices such that Tr p = 1 and Tr pCa = Va, and let 
7^+ be the subspace of V where p is positive-semidefinite. First we note that "P^ is convex: 
Let pi and p2 be elements of V+ and let p(t) = tpi -|- (1 — t)p2,t G [0, 1]. The matrix p{t) 
is clearly Hermitian, properly normalized, and satisfies the constraints. It is also positive 
semidefinite, since for any x, x ■ p(t)x = t(x ■ pix) + (1 — t)(x • P2x) > 0. Thus p{t) is an 
admissible density matrix as desired. 
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Since S* is a strictly concave function defined on a convex space, it cannot have more 
than one stationary point. For if we had two stationary points we could vary p along a line 
between them, and since dPS/dt^ < we could not have dS/dt = at both. 

If there is any p which achieves the maximum value of S", it must be a stationary point 
of S. The only conceivable alternative would be for p to lie on the edge of V+ in "P, but that 
is excluded by the following argument: Suppose p is on the edge of V+ in V. Then there is 
some direction 5p such that p' = p + tSp is positive-semidefinite for arbitrarily small positive 
t but not for negative t arbitrarily close to zero. That can happen only if p has arbitrarily 
small or zero eigenvalues pi which become negative as t becomes negative. For this to happen 
these eigenvalues must increase from arbitrarily small values as t becomes positive. As t — > 
these eigenvalues dominate everything else in dS/dt = J2i (— Inpi), so dS/dt — > +oo. Thus 
for sufficiently small but positive t, S{p') > S{p), so p does not maximize S. 

Whether a maximum-entropy point exists at all depends on the constraints C^- First 
suppose there is a constraint on the energy Tr pH = Eq and no other constraints. In a 
normal system there will be a state with the maximum entropy.§ Now suppose that there 
are additional constraints. They further restrict the space of allowable p and thus can only 
reduce the achievable values of 5, so there is at least a supremum of possible values of S 
in the constrained system. We believe but have not been able to prove that in our problem 
there will be a p which achieves the maximum value of S. In the numerical work described 



in Sec. VII we have always succeeded in finding a solution. 



B. The form of p 

We must find the unique state p which satisfies our constraints and which gives 
Tr^plnp = for any Sp which maintains the constraints. This means that we are con- 
cerned with Sp such that 

TtSp = (2.6a) 
TT6pCa = (2.6b) 

for all a. 

Now we can treat the space of Hermitian operators as an inner product space with 
{A, B) = Tr AB. In this language, Eqs. (^I6|) imply that the possible Sp span the space 
orthogonal to the identity and to each of the Ca- Since In p must be orthogonal to all the 
6p, it must be composed only of / and the Ca, so we can write 

In p = const + ^ faCa (2.7) 

a 

for some coefficients /„. Since Tr p = 1, we can write 



^But note that if the density of states rises exponentially or faster with energy, then the entropy 
can be unbounded. 
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focCcx 

P^—y^. (2.8) 

Our goal is now to determine the coefficients fa so that the constraints are satisfied. 
We can define a grand partition function, 

Q = TreS^"^". (2.9) 

Its derivatives are 

dQ 



so 



^ = TrC„eS^"^«, (2.10) 

(C«) = -^lng. (2.11) 

We have the usual thermodynamic formula for the entropy, 

S = -(Inp) = Ing - ^/,(C,) . (2.12) 
Differentiating this we find 

dfa dfo. 

Now we specialize to the case where one of the constraints is just the Hamiltonian. The 
corresponding coefficient is written —f3, and we have 

p^i^e-^^+S/.c«_ (2.14) 

If we vary the coefficients in such a way that (H) = E changes but the other (Cq) = Va 
remain fixed we see that 

dS = -Y,fad{Ca) = (3dE. (2.15) 

Thus /3 = dS/dE and so the coefficient j3 has the usual interpretation as the inverse tem- 
perature. 



III. ONE SCALAR FIELD, ONE DIMENSION. 



Now we will restrict ourselves to a theory consisting only of gravity and one massless 
scalar field. Presumably such a system has enough richness to contain the usual black hole 
information paradox, so nothing important is lost by making this reduction. We will also 
begin here by working in one dimension. We will put our entire system in a box of length L 
and require that all deviations from the vacuum are in the region from to Lin. Later we 
will let the overall box size L go to infinity while Lin remains fixed. The inside region will be 
[0, Lin] and the outside region will be [Lin, L]. We will use the usual scalar field Hamiltonian, 
which in classical form is 



H 



nix) 




dx . 



(3.1) 
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IV. THE SOLUTION DENSITY MATRIX IS GAUSSIAN 



In our problem, the constraints are the energy bound, 

Tr pH^Eo, (4.1) 

and the "vacuum-bounded" condition, 

TrpO^^iOmO), (4.2) 

where Of is any operator which is constructed out of the fields (f){x) and 7r(x) in the outside 
region. The solution must have the form 

poce-'^^+S^-^r (4.3) 

We now show that fa is nonzero only for those operators which are quadratic in the 
fields. 

Suppose that we wanted to solve a different problem in which we cared only about the 
constraints involving the quadratic operators. We would have the energy bound and the 
constraints 



IVpQ^ = (0|Q«|0), (4.4) 

where Q2 runs only over quadratic operators (f){x)(f){y) and 7r(x)7r(7/). (The operators 
(l){x)7i{y) vanish automatically by symmetry under — > —0.) The solution to this problem 
has the form 

p'cxe-^^+^f^"^^. (4.5) 

Now p' is a Gaussian operator, i.e. {(Pi-)\p'\4''{-)) is a Gaussian functional of the values 
of and 0'. Let p'2 = Tri p'. The trace is just a set of Gaussian integrals, which means 
that the resulting p'2 is also a Gaussian. Because H is quadratic, the vacuum p^ — |0)(0| is 
Gaussian, and so is its trace p^ — Tri p^^. Now by construction we have Tr p2(f){x)(f){y) — 
Tr P2^'^(f){x)(f){y) and Tr p'yn{x)7r{y) = Tr p™^7r(a;)7r(|/). These conditions are sufficient to fix 
the coefficients in the Gaussian p'2, and thus to show that in fact p'2 and p^'^ are the same 
Gaussian, i.e. that Trip' = Trip™*^. 

Thus p' satisfies all the constraints of the original problem. Since only one p can have 
these properties it follows that p — p' and thus that Gaussian solution p' is the correct 
solution to the original problem. 



V. THE DISCRETE CASE 



We now approximate the continuum by a one-dimensional lattice of coupled oscillators, 
with a classical Hamiltonian 

//=^(P,-P, + x.Xx) . (5.1) 
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The simple kinetic term in Eq. ( |5.1| ) corresponds to choosing oscillators of unit mass, re- 
gardless of how densely they are packed. In terms of scalar-field variables this means that 
Xi = \/a(l)i and Pi = 7ii/y/a where a = L/{N + 1) is the lattice spacing, 0j is the average of 
(f){x) over an interval of length a, and VTj is the total momentum 7r{x) in the interval. 

The matrix K gives the couplings between the oscillators and represents the d(j)/dx term 
in the scalar field Hamiltonian. To approximate the continuum with the zero-field boundary 
condition we will imagine that we have oscillators located at the points 1/A^+l . . . A^/A^-|-l. 
and the end oscillators are coupled to fixed-zero oscillators at and 1. Then 



K= Q -g ■ ■ ■ (5.2) 



where g = {N + VfjL^. 

We are trying to satisfy the constraints 

TrpH = Eo (5.3a) 

Tt pXiXj = {0\xiXj\0) (5.3b) 

Ti pPiPj = {0\PiPj\0) , (5.3c) 

where i and j run over the oscillators which represent the outside region. We have shown 
that the solution will have the form 



/2g 


-9 


■ \ 


-9 


2g -g 







-9 ■ 












V ■ 




■ / 



p oc e 

We can write this in a more familiar form as 

p oc e-^^' , 

where H' is a fictitious Hamiltonian for these oscillators, 

1 



H'=- 
2 



■ A^'x) 



with 



and 



r 



K' 





' ] 












K'22 



(5.4) 



(5.5) 



(5.6) 



(5.7) 



(5.8) 



Here Ku and K12 are the sections of the original coupling matrix and T22 -^22 
and (3 can be adjusted in an attempt to meet the necessary conditions. This gives us one 
number, /3, and two symmetric Ai'out-by-A^out matrices, and T22, that we can adjust. The 
constraints involve one scalar constraint, for and two iVout-by-Ai"out symmetric matrices 
of constraints, for XiXj and PiPj- There are equal numbers of equations to satisfy and free 
parameters to adjust, and so, if we are lucky, we will be able to find a solution. If we do 
find a solution, we know it is unique from the arguments of Sec. [II A, 
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A. Computing the expectation values 

To actually solve these equations we will need to compute the expectation values of XiXj 
and PiPj given the density matrix p oc exp{—(3H'}. We can compute them in the usual way 
as derivatives of the partition function 

Q = Tre~^'K^M'^^'^^''+^M''^'^''^'') , (5.9) 



so 



To compute the partition function we will first find the normal modes of the fictitious 
classical problem with this Hamiltonian, and then treat these as independent oscillators 
which we will quantize. The Hamiltonian H' gives rise to the equations of motion 

(5.11) 

so we look for eigenvectors x" that satisfy 

T'K'x" = a;^x" . (5.12) 

The eigenvectors will be complete, so that we can define new coordinates Za via 'X. — J2 ^ax", 
which will then obey the equations of motion 

-^ = -u;lza. (5.13) 

We can choose the norms of the eigenvectors so that x" • X'x^ = uj'^5af3 and group the 
eigenvectors into a matrix V via V^q = x'^. Then we will find that = or 

K' = V-^^n'^V-^ and T' = VV'^ (5.14) 
where = ^^a^a/s- We can then substitute 

x = 1/z P^ = V-^^^P^ (5.15) 



into H' to get 



i/' = ^ (p, • p, + z • n'^z) 

-lUPa+^yj^JlH'a- (5.16) 
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This is the Hamiltonian for a set of (fictitious) uncoupled oscillators with frequencies cjj, 
and thus Q = Tie^^^ will be the product of partition functions Qa = Tre^*^^" for the 
individual oscillators. The partition function of a single oscillator is easily found to be 

and thus 

In Q = ^ In Q« = - 5^ In sinh j . (5.18) 



Differentiating gives 



dua 2 2 

Letting A^^ = <^q, this becomes 



InQ = — — coth — a;^ • (5.19) 



^ InQ = --^coth^a;^- (5.20) 



dAa Aoja 2 



We can compute the effect of a change to K[- or T^^ on the eigenvalues, 



ViaVja (5.21) 



and 



"^"^^ V-'V-'Aa. (5.22) 



Finally, 



(XiXj) = Yl ^^iaVja coth (5.23a) 



a 



(^^^.) = E V^.^ coth ^o;. . (5.23b) 

These expressions must have the same values when V and the u are the eigenvectors and 
frequencies of T'K' as they have in the vacuum, i.e. when V and cu are the normal modes 
and frequencies of the vacuum and (5 — oo. 



VI. THE FORM OF THE NORMAL MODES 

We are trying to find the normal mode eigenvectors x° and frequencies uJa that will 
satisfy our constraints. If we take the uia as given, we can learn about the as follows: 
The vectors x"^ satisfy the equation 
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(6.1) 



with 



T' 











^22 , 



(6.2) 



and 



K' 



Ku 



Vo 



-9 



K' 



0\ 



22 



(6.3) 



Let us introduce the convention that Latin letters from the start of the alphabet, a,b,c, . . . 
range over only "inside" oscillator indices, Latin letters from the middle of the alphabet, 
i, j, k,. . . range over only "outside" oscillator indices, and Greek letters range over all indices. 
Writing the eigenvalue equation out in components. 



T' K' r" 



Taking only the inside components of the eigenvalue equation, we see that 



That is to say. 



(6.4) 



(6.5) 



-9x1 + 2gx^ - gx^ 



^a^2 



(6.6a) 
(6.6b) 



(6.6c) 



Taking a;^ fixed, there are N^^ equations involving Ni^ + 1 unknown components of x". 
However, the equations are invariant under a uniform rescaling of x*^. Thus these equations 
fix for = 1 . . . A^in + 1, except for normalization. The equations are readily solved, and 
the solution is 



< = sin Ma, 



(6.7) 



where 



cos k' 



1 _ ^ 
2^' 



(6.8) 



and N'^ is an unknown normalization factor. Here k'^ and N'" can be complex, but x" must 
be real. When k'^ is real we will write ka — k'^ and — N'^ and call this a "normal" mode. 
When k'^ is complex we can write 
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k'^ = n + ika (6.9a) 
K = (6.9b) 
-)^'~^N^smhfiko, (6.9c) 



with ka and Na real. We will refer to these as "abnormal" modes. 

As before, we group the eigenvalues into a matrix, V^a = and arrange the normaliza- 
tions so that K' = V-^ AV'^ and T' = VV^ . We know the value of T' if one of its indices 
is an inside index, 



Jab — J-ab 



VaaVba (6.10) 



and 



Q = Taj=VaaV,^. (6.11) 

From Eqs. (|6.7|) and ( |6.10|) we get 

^ab = E sin k'^a sin k'J . (6. 12) 



Since Vat.^+i^o = N' smk'^{Nia + 1), we can use Eq. ( |6.11|) to extend the range of Eq. (|6.12|) 



to 1 . . . A'in + 1 except that the equation does not hold when both a and b have this value. 

A similar calculation can be done for V~^, the inverse of the eigenvalue matrix V. In 
this case we will find that 

V-,' = Ksmk'^a = Vaa. (6.13) 

This has the same form as the equation for V, but applies only for a = 1 . . . Ni^. That 
makes V~^ less useful than V for establishing a connection between the inside and the 
outside region, and we will not use it further. 



VII. NUMERICAL STUDIES 



We have solved numerically the set of nonlinear equations, Eqs. ( p.23| ), with (3 fixed and 



various values of Lin and A^. The problem is one of solving Aout(Aout + 1) simultaneous 
nonlinear equations for Aout(Aout + 1) parameters. In general such problems are quite 
difficult to solve, even if we know that there is a unique solution. Here we used the Powell 
hybrid method pH|. If there are no local minima of the rms error in the function values. 



this method converges from any starting point. Fortunately this appears to be the case 
in our problem. However, Powell's method often converges quite slowly for large systems, 
requiring many thousands of iterations to make progress. This has limited our numerical 
solutions to problems with no more than about 30 oscillators. The codes were written in Lisp 
and executed on DEC Alpha workstations. The calculations were done using an arbitrary- 
precision floating point package. The results presented here were computed using at least 
38 digits of precision and the solutions found were accurate to at least 21 significant digits. 

To understand the numerical solution we look at the normal mode frequencies and the 
forms of the normal modes. When the mode is "normal" (i.e. real k'^ in Eq. ( |6.7| )) the mode 
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Mode 1 Freq 3.2798 ^ 
Mode 2 Freq 6.8002 -+- 
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FIG. 1. Modes and frequencies for the "normal" modes of a system with L = 1.0, Lin = 0.5, 

13 = 2, N = 6. 
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FIG. 3. The second normal mode, rescaled to the continuum normahzation, for various different 
number of oscihators, with L = 1.0, Lm = 0.5, j3 = 0.5. 

is a sine wave in the inside region. When the mode is "abnormal" it is essentially a growing 
exponential. Typical modes for a small number of oscillators are shown in Figs. |l| and |^. 

As N becomes large, each "normal" mode and its frequency smoothly approach a limit, 
providing that we use a normalization appropriate for the continuum, which means that 



each mode must be rescaled by \J{N + i)/L. See Fig. |^. As N increases, each abnormal 
mode and its frequency undergo a smooth evolution, until at some point it disappears from 
the set of abnormal modes and is replaced by a normal mode with very similar form in the 
outside region, as shown in Fig. ^. Because of this behavior, we believe that if we could 
solve the continuum behavior directly we would find just the "normal" modes. 

The most striking result of the numerical solutions is that the wavenumbers of the "nor- 
mal" modes are quite evenly spaced. The larger the energy, the more accurate is this 
approximation. In Fig. |^, we show the wavenumbers for L = 1.0, Lin = 0.5, (3 = 0.5 which 
gives an energy of about 0.636. Even at this low energy the fit is good to within a few 
percent of the typical wavenumber. For larger energies the points will lie correspondingly 
closer to the line. 

In the L — s> oo limit, our problem has only two dimensionful parameters, Lin and L'o- 
Thus there is only one dimensionless parameter, LinL^o, that characterizes the problem. 
In the 3-dimensional black hole problem the equivalent parameter is REq, which for the 
parameters in Eq. ( |1.3| ) is about 10^^. Thus for application to black holes we are interested 
only in very large values of L^^Eq for which the linear approximation for the "normal" modes 
will be very good. 

Since the "normal" mode wavenumbers extend up to ~ vr, the number of "normal" 
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4 osc, ev#4, fteq 65.8392 
6 osc, ev#5,/freiq 46.8216 -+-- 
8 osc, ev#6/ freb 40.0010 -Q-- 
10 osc, ev#7, fred 36.7291 x 
12 osc, ev#e, freq\35.0883 - 
14 osc, ev#9, freq 34.4434 -j^ 
16 osc, evflO, freq 34.7680 o 
18 osc, ev#11, freq 37.6197 ■ 




0.00 
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0.50 



0.75 
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FIG. 4. An abnormal mode for L = 1.0, Ljn = 0.5, /3 = 0.5. As A'^ is increased the mode 
changes smoothly until at AT = 18 there is no corresponding abnormal mode, but instead one of 
the new normal modes has a very similar form but different normalization. 



2.5 



^ 1.5 



0.5 



numerical result o 
best-fit line — 



1 2 3 4 5 6 7 

Mode number 



10 11 12 



FIG. 5. The numerically-computed wavenumbers compared with the best-fit line through the 
origin for L = 1.0, Ljn = 0.5, /? = 0.5, AT = 20. 
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modes will be given an integer A^norm ~ vr/fci where ki is the interval between wavenumbers. 
In addition there will be — A^norm "abnormal" modes, with frequencies u > 2{N + 1)/L. 
For N ^ L/f3 these modes do not contribute to the entropy, because they are exponentially 
suppressed. 

VIII. THE FIRST OUTSIDE OSCILLATOR 

In an attempt to minimize our reliance on numerical results we will use only the statement 
that the wavenumbers are evenly spaced. This leaves one free parameter, the spacing ki, 
which depends on the energy Eq. To fix ki we will look at the correlator of the position of 
the first outside oscillator with itself, which from Eqs. ( ^.23a| ) and is 



(x7Vi„+iXjVi„+i) = E ^coth (^^^a^ (^^L+i)' 



= E ^ coth (^^u^^ N'J sin^ k'^{N,^ + 1) . (8.1) 

The important point here is that we know for fi up to A/'in + 1, and we know {x^Xy) for /i 
and V down to A^in + 1. Thus taking n = u = A^in + 1 gives the unique correlator for which 
we know the components that go into the expression for the correlator while also knowing 
that the correlator must have the same value as in the vacuum. The same argument does 
not work for (P/v.^^+iP/v.^^+i), because we know only up to = A^in and not /i = A^in + 1. 

By computing this correlator in the vacuum and in a vacuum-bounded state for a given 
f3 we can fix the spacing of the normal mode frequencies. In a high energy state, (3 will be 
small and so there will be frequencies with 1- For such a frequency, 

coth ~ -T— > 1 • (8.2) 

In order to keep (x^.^^^) the same as in the vacuum we must require 

sin2A;^(Arin + l) < 1 (8.3) 

to cancel the large coth, and consequently 

ka,{N;r, + 1) ^ nvr (8.4) 

for some integer n. This means that ka is close to one of the wavenumbers appropriate to 
the problem with a rigid boundary at Li^. 

A. The vacuum 

First we will compute the correlator for the vacuum, (x^.^ _,_]^)o. Here /5 = oo so the coth 
term drops out, giving us 
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(^k+i)o = E sin' /Cn(A^in + 1) . (8.5) 



For the ground state of N oscillators we have 



and 



which give 



iV+l . K 

= 2 ^ sm y (8.7) 



(^Afin + l)0 



^ JV+1 

2(iV+l)„^i(A^+l)sin^- 



We would like to evaluate this expression in the N ^ oo hmit, with Niji/N fixed. There 
is a prefactor L/{N + 1), which goes to zero in this limit, but that is just an artifact of 
the conventions we have used for the discrete problem, and will appear in the finite-energy 
vacuum-bounded states as well. 

We will expand the sum for large N. There is a logarithmic divergence. We are interested 
in the In N term, and in the constant term, but we will ignore any terms of order 1/N or 
below. Accordingly we will take N + 1 a,s N wherever it occurs. 

First we expand the numerator using sin^x = (1 — cos 2a;) /2. The first term gives 

1 ^ 1 

- y ■ (8.9) 

2 ^ TV sin — 

^ n=l 2N 

We would like to turn this sum into an integral in the N ^ oo limit. However, we must 
first subtract out the divergent part. Accordingly, we write this as 

1 1 1 / 1 2 \ , , 



The first term gives 



-(In AT + 7) (8.11) 

TT 



where 7 is Euler's constant. The second term is finite and can be converted to an integral 
in the A^ — > 00 limit, to give 



1 /■! / 1 2 \ 



2 lo 



dx — = - 



sin ^ TTX J 



1 





1 , 4 



1 






TTXy 




In 




tan — 


TT 




KX 


4 /. 



(8.12) 



= -ln-. (8.13) 

TT TT 



The remaining term in the sum is 
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^ n=l 2N 

To compute this we use l/sinx = cscx = 1/x + x/6 + 7x^/360 + ■ • ■ to get 

1 ^ 2n7riVi, / 2 vr n 77r n3 \ 

— > cos H 7 + ■ ■ ■ . 8.15 

N \7m 12 2880 J ^ ' 

The first term can be summed, 

TT^ n TT \ N J ^ ^ 

The rest of the terms do not contribute. Because of the oscillations of the cosine, 
X^i^ '^'^ cos(n7rA^in/A^) goes as A^'^ rather than A^^+^ and thus is killed by the correspond- 
ing N^~^^ in the denominator. 

Putting all the terms together we find that 



In + 7 In — sin 



0(1 



71 N 



(8.17) 



B. Using the even-spacing ansatz 

Now we will compute the same correlator in the vacuum-bounded system, using the 
ansatz that the wavenumbers are multiples of some spacing ki, which gives A^norm ~ vr/fci 
normal wavenumbers. This is equivalent to saying that the normal wavenumbers are those 
appropriate to a problem with a rigid boundary at distance Ljn' = {ji /ki){L/{N + 1)) ^ 
L(A^norm/^)- There will also be a contribution from the abnormal modes. The contribution 
from each mode is positive, so by taking only normal modes we will find a lower bound on 
{xlf. j^i). Using this bound we can derive an upper bound on L-^^ and thus on S{E). 

We will again work with A^ oscillators for A^ large. We will ignore all terms of 0(1/A^) 
or 0(l/A^norm)- Accordingly we will take Lin'/^norm = L/N and 

k^ = nki = = . 8.18 

N N + 1 

^ ' norm ^ ' norm ' ^ 



As before we have 



2Ar . /c„ 2Ar . vrn 
ujn = —j— sm -— = — — sm — — . (8.19) 

L Z h ZiVnorm 



Now we can derive the normalizations of the normal modes: We know from Eqs. (|6.10|) 
and (|6.11|) that if a is an inside oscillator, then for any /i. 



Y.^a.V^u = 6a^. (8.20) 
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so 

Nn orm 

Safi = X! sin kna sin A;„/i + abnormal modes. (8.21) 

n=l 

Later we will find that A^norm is close to A^in- Using this, and since we are working in the 
limit where A^in is large, we can choose a such that 

1< a < 2Ari„ - N,,orm ~ A^i^ . (8.22) 

Each abnormal mode a contributes 

A"^ sinh katt sinh k^fi = (5abn (8.23) 

to the right-hand side of Eq. ( ^.211 ). However, this same mode contributes 



L ^^, sinhH,(Ari, + l) 

AN cosh(A;„/2) ^ 
to the sum for {x%.^^^) in Eq. (|]l|). Since (a;^.^^+i) = L/(27rA^)[ln A^ + 0(1)] it follows that 

^^^ sinh^ A;,(A^in + l) < InAT 

cosh(A;j2) ~ 27r ^^'^^^ 

and thus that 

sinh /c«a sinh cosh(fca/2) , 

Oabn < o ^ ^ llliV . (8. 2d I 

27rsinh2A;«(Ari, + l) ^ ^ 

The exact values of the k^ for the abnormal modes vary with A^. For successive A^ 
values, each k^ decreases toward until the corresponding mode converts to a normal mode 
as described in Sec. |VI1| . By avoiding the points where these "conversions" are about to 
take place it is possible to find a sequence of values for A^ which have all fc^ ~ 1- For such 
values. 



From Eq. (8.22) the exponent is ^ 1, so 5abn is exponentially small, and the contribution of 



the abnormal modes to (|8.21|) is suppressed 



Thus we ignore the abnormal modes in Eq. ( p.21[ ), multiply by sinfcm/i and sum over n 
to get 

Nn orm 

sin kmCt = ^ ^ sin kmf^ NI sin k nOj sin /Cfi/x 

/i=l n 

-^norm 



EAr2 • ; • vrmu . Tina 
N„ sm kntt > sm -^j — ^ sm ^ 



EiV^infc„a- ^— + (8.28) 



A^ + 1 

j2 norm \ ^ 



K : sinfc^a, (8.29) 
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from which we conclude that 



+ 1 iV„ 



(8.30) 



Thus we have 



sin k„Njr 



n=l 



coth 



f8.311 



norm^n 



Now define A by 1 — A = A^m/iVnorm ~ -^in/-^^in', SO that sin^ fc„A^in = sin^n7r(l — A) 
sin^ nvrA. Then put in the definition of Un from Eq. ( ^.191 ) to get 



)> — Y 



sin^ nvr A 



sm 



2iVn 



, ( pN nvr 
X coth —— sm 



2N^ 



(8.32) 



In the large limit, the argument of coth becomes PttuN / {2LNj^orra) ~ P^^n/ {2LiJ) = 
7™/(2r') where r' = L^^'T = LiJ / (3. Now if it were not for the coth term, we would have the 
same sum as in Eq. ( ^.8|) with iVnorm for except in the prefactor. Thus we will exchange 
cothx for 1 + (coth X - 1) = 1 + 2/(6^^ - 1) in Eq. to get 



> 



L 



2nN 



In A^norm + 7 + lu ( - SlU TT A 



where 



B 



TT 



sin tittA 



sin p 



vrn/r' 



(8.33) 



(8.34) 



Now we again use 1/ sinx = 1/x + x/6 + ■ ■ ■. The first term has no A^norm dependence 
and we can extend the sum to oo. In the next term, the sum is cut off by the exponential in 
the denominator, leading to a term of order (T'/A^norm)^- Further terms have higher powers 
of t'/N^oj.^. In the limit A^norm ^ cc we ignore all these terms, which leaves 



5 = 4E 



sin^ mi A 



\ n [e^ri/r' _iy 



(8.35) 



We are interested in the high-energy limit, for which r' ^ 1. Later we will see that A is of 
order In r'/r' <^ 1. Thus the summand in Eq. (|8.35|) is slowly varying and we can convert 
the sum into an integral. 



5^4 



°° sin TcAxdx 
X (e'^^/^' - 1) 



(8.36) 
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The error in Eq. ( p.36|) is approximately the term that we would have for n = in Eq. 
( ^.35| ). Taking the n — > limit we find that this term has order A^r' ~ (In r')^/r' ^ 1, so 
our approximation is good. The integral in Eq. ( p.36| ) can be done and the result is 

_ -47rr'A 

5 = 27rrA + ln — . (8.37) 

47rr'A ^ ' 



Now we set {x%^_^-^) = {x%.^^+^)o to get 

In A^norm + 7 + 1^ (^^ siu TT + B 



or 



/8 ttN- \ 

<lnA^ + 7 + ln -sin— ^ (8.38) 
Vtt a* / 

, / N sin^\ , / L sin^\ 
5 < In = In ^ . (8.39) 

V^Vnorm Sm TT A / \J^m Sm TT A / 

Note that we have taken the limit N —>■ oo and there is no longer any dependence on A^. If 
we now let L — cxD with Ljn fixed, we can approximate sin(7rLin/L) = tiLi^/L to get 

5 < In .^^^ , = In ~ In ^ (8.40) 

Lin sinyrA sinvrA ttA 

since A is small. 

Using B from Eq. ( |8.37| ) we get 

1 — e"^'^^'^ 1 

27rr'A + ln < In — . (8.41) 

47rr'A " vrA ^ ' 

Thus A < Amax where 

1 At' 

A„,ax = TT^ In ... . (8.42) 

If instead of r' = L^-^T we use r = LjuT we will make an error of order A^ax; which we 
expect to be small. So we ignore the second order contribution and take r' as r in Eq. ( ^.421) , 

^max = In -. . (8.43) 

If we ignore e~^'^'^^™=^ in the denominator we get 

A„,ax = 77^1n4r. (8.44) 

Using this we find that e"^'^^^™'"' = (4r)^^ ^ 1 since r ^ 1, which justifies ignoring this 
term. We will also ignore In 4 by comparison with Inr. Thus we conclude 

^ - 2^ ^ + K^) - 2^ + ^ (Z^) ^'-''^ 



Ur! <U^ + ^ In Li„T + 0(l-] . (8.46) 



and 

2-kT VT. 

The equivalent system is larger by at most a thermal wavelength times a logarithmic factor 
depending on the inside size. 
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IX. PROPAGATION OF BOUNDS 



In the previous section we derived an expression that gives the frequencies, and thus the 
entropy, for a vacuum-bounded system at a given temperature T = 1/(3. Given such an 
expression, we would hke to compute the entropy as a function of energy. But the energy 
is not simple to compute from the frequencies alone.0 However, we can easily compare the 
entropy of the vacuum-bounded system to that of a system with a rigid boundary at Lin 
and the same temperature. To make this comparison at fixed energy instead, we proceed as 
follows. 

Consider the free energy F = E — TS which has dF = —SdT. Integrating gives 

E-TS = - r S(T')dT' . (9.1) 
Jo 

Let 5''''^, T"^^ denote the entropy and temperature of the system with a rigid boundary at 
Lin, and S, T, and so on denote those for the vacuum bounded system. For any quantity 
A let 6A{T) denote the difference between vacuum-bounded and rigid-bounded systems at 
fixed temperature, 6A{T) = A{T) - A'^{T), and 6A{E) denote the same difference at fixed 
energy, SA{E) = A{E) — A^^{E). With E fixed we compare the differences (to first order) 
in the two sides of Eq. (|9.1| ) between the vacuum-bounded and rigid-bounded systems, 

-T6SiE) - ST{E)S = -5T{E)S - 5S{T')dT' , (9.2) 

Jo 

where the first term on the right-hand side comes from the change in the integration limit. 
Thus 

3S{E) = ^ r 6S{T')dT'. (9.3) 



TJo 



X. DERIVING THE BOUND 



Now we apply these results to the case from Sec. [VIII B| where 



with 



and Lin' Lin(l + A) (10.1) 



ALin< -^InLinT. (10.2) 

ZTTJ 



'^It can be done, but since E needs to be renormalized against the ground-state energy of the 
entire system the result depends sensitively on the frequencies and normalizations even for very 
high-energy modes. 
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This says that at any given temperature, the vacuum-bounded system has the entropy 
S{T) of a system of length L[J. Now in a one-dimensional system the entropy density is 
proportional to the temperature, 



S^^ = ^Li„T, (10.3) 
and thus the entropy difference between vacuum-bounded and rigid-bounded systems is 

6S{T) = ^ ALi,T < i In h^T . (10.4) 

6 O 



Using the results from Sec. 0, 

SS{E) < i r 1 In L,XdT' = i (In L,^T - 1) . (10.5) 
I Jo b 

Since we are ignoring terms of order 1 by comparison with those of order In T, we can write 
this§ 

55(E) < ^ In LinT. (10.6) 



So, we conclude that the vacuum-bounded condition closely approximates the rigid box 
of length Lin. For the same energy, the vacuum-bounded condition allows slightly more 
entropy. The entropy difference grows at most logarithmically with the number of thermal 
wavelengths in the inside region. Since the total entropy at a given temperature is linear in 
T, the ratio of dS to S goes to zero at high temperatures. 



XI. DISCUSSION 

We have introduced a new way of specifying that matter and energy are confined to a 
particular region of space. Rather than giving a boundary condition per se, we specify a 
condition on a density matrix describing the state of the overall system. We require that any 
measurement which does not look into the inside region cannot distinguish our system from 
the vacuum. This avoids certain difficulties such as the Casimir energy that results from 



the introduction of a boundary and the geometric entropy p3|,^ that results from ignoring 
part of a system. For these "vacuum-bounded" states, we consider the problem of finding 
the maximum-entropy state for a given total energy. This is analogous to the problem of 
finding the thermal state in a system with a rigid boundary. 

Unfortunately, the vacuum-bounded problem is more difficult than the analogous prob- 
lem with a rigid boundary and we must resort to working in one dimension and to numerical 
solution on a lattice. From the numerical solution we justify the ansatz that the continuum 



^^It has happened that 6S{E) and 6S{T) are approximately the same, but that is a particular 
property of the system at hand. For example, if A were a constant we would have 6S{E) = 
1/T cALi^T'dT' = 2cALinT = 26S{T). 
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wavenumbers are evenly spaced in this problem. Using this ansatz we compute an upper 
bound on the entropy of a vacuum-bounded state, and show that for high energies {ER ^ 1) 
the entropy approaches that of a system with rigid boundaries. Of course this is what one 
would expect for a system whose typical wavelengths are much shorter than the size of the 
inside region. 

To apply this result to an evaporating black hole we look at the state produced by the 
black hole after evaporation . Since our calculation was one-dimensional we must assume 
that the similarity between the vacuum-bounded state and the thermal state with a rigid 
boundary extends to 3 dimensions. Then we infer that very little entropy can be emitted 



in the final explosion, confirming the results of Aharonov, Casher and Nussinov |T^ and 
Preskill [H. For example, a black-hole formed in the big bang with mass of order lO^^^f 
would be evaporating today. During its life it would have radiated entropy S ~ 10^^. Now 
we assume that the entropy of the final explosion has energy E ~ lO^^erg contained in radius 
R ~ 10~^^cm as in Sec. [1 B[ and that the maximum entropy is not too different from that 



of a spherical box, in accord with our 1-dimensional result. Then we find that the final 
explosion can emit only entropy S ~ 10^'', which is a factor of 10^^ too little to produce a 
pure state. 

This argument means that either a black hole must not evaporate completely but rather 
leave a remnant or remnants, that information must be lost, or else that the Hawking 



radiation is not exactly thermal, even at very early times [12 
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APPENDIX A: CONCAVITY OF THE ENTROPY 

Lieb [23] showed that 5* is always (downward) concave. Here we obtain the same result 



by a different technique and show also that the concavity is strict, i.e. that (PS/dt^ never 
vanishes. 

We consider S = S{p') = —Trp' In p' with p' = p + t6p for some 6p with Tr 6p = 0. Then 

— = -Tropin p' (Al) 
at 



^^These routines and many other are now available on a Numerical Recipes CDROM [27| 
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and 



Tr Sp— In p' . 



df^ dt 
To expand this we use the formula 
d 



, In 

dt 



A = ds (/ - s{A - I))-' ^ (/ - s{A - I))-' 



which can easily be derived as follows: Let B = A — I. We expand 



fi2 53 

In A = ln(/ + B) = B + — + — + - 



and differentiate to get 



d ^ ^ dBl dB „ dB\ 
— In A = — — —B + B— 
dt dt 2 \ dt dt J 

1 fdB „ dB ^dB\ 



+ - —B' + B—B + 5^— + 
3\dt dt dt J 



We observe that this is 



1 , dB fdB dB\ 2 fdB . dB .dB\ 

ds— + s —B + B— ]+sU —B^ + B—B + B^ — 
at \ at at J y at at at J 

'l + sB + s^B' 



£ ds (l + sB + s^B^ + ■ 

C ds {I - sB)-' ^ {I - sB)-' 
Jo dt 



as desired. 

Using Eq. ( |A3| ) we find 

d^__ 
dt^ ~ 



C dsii 5p{i - s(p - i)Y^5p (I - s(p - nr^ 

Jo 



^iXisf 

Jo 
where 

xis)^i5py/'ii-sip-i))-\6py/\ 

Since X{s) is Hermitian, TrX(s)^ > with equality obtained only for X{s) 

dt^ - 

with equality only if X(0) = for all s. But X(0) = Sp, and thus 



dt^ 

for any nonzero 6p. 
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